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ABSTRACT 


The National Aeronautics and Space Administration Armstrong Flight Research Center (Edwards, 
California) Flight Loads Laboratory is a unique national laboratory that supports thermal, mechanical, 
thermal/mechanical, and structural dynamics research and testing. A Proportional Integral thermal control 
system was designed and implemented to support thermal tests. A thermal control algorithm supporting a 
quartz lamp heater was developed based on the Proportional Integral control concept and a linearized heating 
process. The thermal control equations were derived and expressed in terms of power levels, integral gain, 
proportional gain, and differences between thermal setpoints and skin temperatures. Besides the derived 
equations, user’s predefined thermal test information generated in the form of thermal maps was used to 
implement the thermal control system capabilities. Graphite heater closed-loop thermal control and graphite 
heater open-loop power level were added later to fulfill the demand for higher temperature tests. Verification 
and validation tests were performed to ensure that the thermal control system requirements were achieved. 
This thermal control system has successfully supported many milestone thermal and thermal/mechanical tests 
for almost a decade with temperatures ranging from 50 °F to 3000 °F and temperature rise rates from -10 °F/s 
to 70 °F/s for a variety of test articles having unique thermal profiles and test setups. 


NOMENCLATURE 


A/D 

analog-to-digital converter 

AGP 

acquisition and control processor 

acoeff 

A coefficient 

ANSI 

American National Standards Institute 

AFRC 

Armstrong Flight Research Center 

bcoeff 

B coefficient 

c, C 

controlled variable, skin temperature, °F 

C/C 

carbon/carbon 

CLM 

closed-loop map 

CMC 

ceramic matrix composite 

C/SiC 

carbon silicon carbide 

d 

number of delay cycles 

DACS 

Data Acquisition and Control System 

D/A 

digital-to-analog converter 

e, E 

thermal setpoint minus skin temperature 

F 

number of cycles per second 

FLL 

Flight Loads Laboratory 

GH 

graphite heater 

I 

current signal 

Im (z) 

imaginary axis in the z-domain 

K, 

integral gain coefficient 

K P 

proportional gain coefficient 

n-1, n, n+1 

discrete time 

NASA 

National Aeronautics and Space Administration 

PCC 

power control cabinet 

PCM 

profile control map 

PI 

Proportional and Integral 

QL 

quartz lamp 

QLGH 

quartz lamp and graphite heater 

Specld 

specimen identification 

Re (z) 

real axis in the z-domain 
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SCR 

silicon controlled rectifiers 

t 

time 

TC 

thermocouple 

TCS 

thermal control system 

THRM PET 

thermal profile elapsed time 

r, R 

reference input, thermal setpoint, °F 

u, U 

controller output, power level 

Z 

z-transform variable 

a 

integral sensitivity factor 

P 

proportional sensitivity factor 

At 

time interval, second 


INTRODUCTION 

A Proportional Integral (PI) closed-loop thermal control system (TCS) was built in-house to support 
large-scale heating tests at the National Aeronautics and Space Administration (NASA) Armstrong Flight 
Research Center (AFRC) Flight Loads Laboratory (FLL) in Edwards, California. A complete TCS includes 
mathematical models, algorithm, hardware, and software. This report emphasizes the mathematical models 
and algorithm; it does not describe the TCS supporting hardware and software. The TCS consists of two basic 
parts: controller and heating process. The controller includes a real-time operating system acquisition and 
control processor (ACP), input modules or analog-to-digital (A/D) converters, and output modules or 
digital-to-analog (D/A) converters. The heating process contains a power control cabinet (PCC), quartz lamp 
(QL) heaters, and control thermocouples (TC)s. 

The design of the TCS involved a popular PI control model and a linearized heating process model (ref. 1). 
The equations that represent the models were mathematically derived at first in continuous-time domain, and 
later changed to discrete-time domain orz-domain (ref. 2) to be implemented in a real-time operating computer 
system. Then, the two models were combined to derive the closed-loop TCS transfer function and stability 
region. With some reasonable assumptions to simplify the mathematical derivation process, a closed-loop 
TCS transfer function in the z-domain was achieved, and a stability region was determined. A method to 
calculate the integral gain K t and the proportional gain K p was developed so that the gains were within the 

stability region. The derived thermal control equations and thermal maps that contain important information 
were used to implement and to add new capabilities. The thermal control algorithm was written in American 
National Standards Institute (ANSI) C programming language to build the TCS that met all of the QL thermal 
control system requirements (ref. 3). The executable code of the thermal control algorithm was loaded into 
the ACP (ref. 4), a component of the FLL Data Acquisition and Control System (DACS). Since the PCC is 
the main component of the heating process, the operation of the PCC zero crossover is described to explain 
the role of the PCC in the derivation of the thermal control equations and its function of firing the QL to 
provide radiant heat. 

In 2005, graphite heater (GH) power control cabinets were integrated into the TCS to support heating tests 
with temperatures up to 3000 °F. In 2006, GH closed-loop thermal control capability was added in the TCS 
using the derived thermal control equations and software control logic based on test data and the GH high 
thermal inertia characteristics. The GH closed-loop thermal control software was tested in a small test chamber 
that can support up to 2 zones. When finished, the GH closed-loop thermal control had all the capabilities 
applied to the QL closed-loop thermal control. 

In 2008, GH open-loop power level capability was added to allow the TCS to fire the GH with predefined 
power levels programmed in the thermal maps. Since GH has high thermal inertia, control TC was required 
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in each zone to monitor the highest temperature limit for safety. The GH open-loop power level was also 
achieved using the small test chamber. 

In 2010, both QL and QH were installed in a test setup for very high temperature tests. The TCS produced 
high power levels to provide heating conditions with temperatures up to 3000 °F and with high temperature 
rise rates. 

From 2003 to 2010, this TCS had performed many large-scale heating tests with either QL, or GH, or 
both QL and GH (QL GH). The TCS had successfully accomplished the criteria of all thermal-structural test 
programs operating in the FLL. With this robust and diversified TCS, the FLL was able to support hot 
structural research for advanced aerospace programs including milestone thermal and thermal/mechanical 
tests that had never been done before. 

QUARTZ LAMP THERMAL CONTROL SYSTEM DESIGN 

The design of the QL TCS involves a controller and a heating process as shown in figure 1 . The controller 
outputs the required power level to the heating process to heat a test specimen following a preplanned thermal 
profile. The algorithm of the closed-loop TCS starts in continuous-time and then converts to discrete-time. 
Since computers run their tasks repeatedly every cycle, the discrete-time model must be developed using the 
z-domain. 


PI Control System Model 

The PI controller was selected in the design of the TCS. This algorithm involves two separate parameters: 
the proportional P depends on the present error, and the integral / depends on the accumulation of past errors. 
The PI controller block diagram is described in figure 2 with the assumption that at the initial time t 0 , 

controller output u(t 0 ) =0. Equation (1) is the mathematical representation of the PI controller; 


lt(t) 


= Ki 



dr + K p e(t ) 


( 1 ) 


in which t is the current time, t 0 is the initial time, T is the variable of integration taking on values from 
time 0 to the current time t, r is the reference input, c is the controlled variable, e = r-c is the difference 
between r and c , u(t) is the controller output at time t, is the integral gain, and K p is the proportional 
gain. By taking derivative of equation (1), the control model can be described by the differential equation (2): 


de 

— = K ie + K p — 


du 

dt 


( 2 ) 


Integrating equation (2) from time t ll _ l to time t n with the time interval dt=t n —t n _ j is shown in equation 

( 3 ): 



du = Ki 



edr + K p 



( 3 ) 
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Equation (3) can be integrated to yield equation (4): 


Ur, 


= £/ n -i + K t [ edr + K p (E n - E^) 

J tn -i 


(4) 


where U n is the controller output at time n, U n _ j is the controller output at time n-1, E n =R tl —C n is the 
difference between reference input and controlled variable at time n, E n _ l = R n _\ —C n -1 is the difference 
between reference input and controlled variable at time n-1, t is the variable of integration taking on values 
from time t n _ l to time t n . \' t " ^ ed T is calculated by linear approximation, and roughly equals the area of the 

surface shown in figure 3, the area of a trapezoid with the bases E n and E n , and the height of At. The 
integral part of equation (4) can be expressed in equation (5): 



edr s At 


E n + En - 1 

2 


(5) 


The following equation is obtained by combining equations (4) and (5). The equivalent discrete-time model 
is expressed below in equation (6): 


At 

U n = U n - 1 + K iY + E n-l) + K p( E n ~ E n- 1 ) 


( 6 ) 


Expand equation (6) and rearrange the terms (see Appendix A), the equation of the discrete-time PI control 
model is derived as equation (7): 


U n = U 0 + KiAt 


i=o 


+K v (E n - E 0 ) 


(7) 


There is always an inherit delay in any control system. Equation (7) allows U n to be calculated from 
U 0 ; therefore, U „ can also be calculated from U n _ d , where d is the number of delay cycles between the 
commanded controller output U n _ d to get the reference input R n _ d at time n-d and its current corresponding 
controlled variable C n at time n. Replacing time n-d for time 0 in equation (7) gives equation (8): 


Un — U n -d + K t At 


I ,-(^) 


j=n-d 


+E v (E n - E n _ d ) 


( 8 ) 


Equation (8) is much more practical than equation (7) since the computer can easily store d values of U 
and E before the current cycle or time n, but not from starting time 0. Let time interval At = 1 s and F be the 
number of cycles per second, then the controller output U n at time n is determined as follows in equation (9): 
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( 9 ) 


r 

n 


lUn-a+K, 

V , (Cn+E n - d \ 

+ K p(E n ~ E n -d ) | 

2-U 1 \ 2 ) 

{ 

j=n-cL 

J 


Equation (9) is the controller output equation derived from the PI control model. Equation (9) is the 
general closed-loop PI control equation and can be used for any PI control system with controller output U. 
For the thermal control system, R is the thermal setpoint, C is the skin temperature, E = R - C is the 
difference between thermal setpoint and skin temperature, and U is the power level. 

Heating Process Model 

A nonlinear continuous-time model for the heating process is described generally as follows in equation 

(10): 


- C 4 ) = /(c,it) (10) 

dt iv 2 

where c is the skin temperature, u is the power level, and f(c,u) is a function of c and u. To the above 
nonlinear system, a simplified linear model of the heating process model described in figure 4 is used and 
expressed mathematically in equation (1 1): 

dc 

— — ac + bu+ q (11) 

dt 

Where a, b, g coefficients are not constants and are functions of outside temperatures. The equivalent discrete 
heating process model is obtained at times n and n+1 with the time interval At = t n+1 — t n as shown in 
equation (12): 


r _ r 
'-•n+l 


At 


aC n + bU n + g 


Move the term C„ +1 in equation (12) to the left side as in equation (13): 


( 12 ) 


C n+1 = (1 + aAt)C n + bAtU n + gAt (13) 

To have equation (13) expressed in a more compact form, the coefficients of equation (13) are defined in 
equations (14), (15), and (16): 


A = 1 + aAt 

(14) 

B = bAt 

(15) 

<1 

II 

(16) 
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Equation (13) can be written as equation (17): 


C n+ i — AC n + BU n + G (17) 

Based upon the inherit delay in any control system, it is assumed that the power level U applied at the 
beginning, n = 0 , does not cause any changes in the specimen skin temperature. 

At n = 0, U 0 = 0 , and C n+l = C n = C 0 , equation (17) can be expressed as equation (18): 


Cq — tICq + G 

Moving all the C 0 terms to one side yields equation (19): 


(18) 


G 


= C 0 


-A C 0 


(19) 


Equation (19) defines a coefficient G based on initial temperature C 0 and coefficient A . Substituting G 
into equation (17) gives equation (20): 

C n+1 = AC n + BU n + C 0 —AC 0 (20) 

Grouping the A terms yields equation (21): 

C n +i = A(C n - C 0 ) + BU n + C Q (21) 

Equation (21) is the heating process equation. 


Closed-Loop Thermal Control Transfer Function 

After deriving the PI control equation (9) and the heating process equation (21), the PI control model and 
the heating process model are combined as shown in figure 5 so that the closed-loop thermal control transfer 
function can be developed. The transfer function provides a basis for determining important control-system 
response characteristics and stability. 

The time interval in the control and the time interval in the heating process are assumed to be equal. 
Rearrange terms in equation (21) (see Appendix B); the relationship between the controlled variable and the 
reference input can be written as follows in equation (22): 


r — r 

'-•71+ 1 '-•n 


(a - BK p +BKi y) + R n ^BKp—BKi y) + BK t At ^ Rj ~ Z Cj + 7 ° 

" " “ \j=o 7=o / 


( 22 ) 
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where A and B represent the coefficients of the test specimen and are functions of temperatures, and 


At At 

lo = Co ~ AC 0 — BKi —Rq + BKi — C 0 — BK p R 0 + BK p C 0 


( 23 ) 


Equation (23) represents a constant at the initial time n - 0 and is defined based on the initial conditions of 
A,B,K t , and K p . 

In order to derive the transfer function of the TCS, equation (22) must be transformed to the z-domain 
(see Appendix C). The equivalent z-transform of equation (22) is described as follows in equation (24): 


C(z) 


1+A-BK, 

= R(z) 

+ C 0 


A t\ , / At\ 

BK i —) + z - 2 [A - BK p + BK t — Jj 

,l BK p + BK i l Ayz- 1 (uK l 


BK ‘ T) 


+ z 1 (/ 0 C 0 ) 


(24) 


where I 0 -C 0 and C 0 are initial condition terms at time n = 0 . To simplify the developing process, these 
terms are assumed negligible and ignored. Equation (24) can be written as equation (25): 


C(z) _ z _1 ( BK p + BKt x) - z" 2 (BKp - BKt x) 

1 - z” 1 (l+A-BK p - BKi x) + z ~ 2 { A ~ BK P + BK i x) 


(25) 


To simplify equation (25), two new values a and fi are introduced as follows in equations (26) and (27): 


40 

<1 

II 

a 

(26) 

P = BK p 

(27) 


where a is called the integral sensitivity factor, and fi is called the proportional sensitivity factor. Both values 
of a and ft must be positive numbers, a > 0.0 and p > 0.0. Considering equation (14), A = 1 + a At is 
dependent on coefficient a, a function of outside temperature. 


The effect of outside temperature is assumed minimal, a = 0 and equation (14) becomes A — 1 . Replacing 
a and [i into equation (25), a simpler form as in equation (28) is achieved: 


C(z) 

R(z) 


1 — z 


- 1 (2-^f) + Z -( 1 -^f) 


(28) 
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Multiplying both the denominator and numerator with z", the transfer function can be written in the form of 
the powers of z as shown in equation (29): 

cm (/»+¥)■>-¥) 

Equation (29) is the closed-loop thermal control transfer function. 

Thermal Control System Stability 

In the z-domain, the control system is stable if the poles of equation (29) lie inside the unit circle in the 
z-plane. Normally, the poles of a control system have two parts: the real part and the imaginary part. This 
thermal control system is designed by assigning specific values to the poles that are on the real axis and on 
the right-hand side inside the unit circle. These poles have the real part and no imaginary part. This assignment 
leads to 0<z 1? z 2 <1 as shown in figure 6, and equations (30) and (31): 


0 < z 1; z 2 = 

( 2 _, 4 ) 2 _ 4(l _, + | )a „ 

Solving equations (30) and (31) (see Appendix D), a TCS stability region is determined as shown in 
figure 7. The selection of the values of a and p in the stability region must be done carefully since these values 
are used in the calculation of K i and K p . If the proportional gain is too high, the system can become unstable. 

In contrast, if the proportional gain is too small the system can become less robust or less sensitive. The 
proportional term of a PI control system contributes the bulk of the output change. The integral term responds 
to accumulated errors from the past and can cause overshooting the setpoint. After tuning the control 
parameters, the values of a = 0.35 and p = 1.1 were selected inside the system stability region. 

Integral Gain and Proportional Gain K 

In order to determine K i and K p gains that are included in the closed-loop PI control equation (9), a 

pre-test open-loop specimen identification (Specld) test is run. The purpose of this test is to compute A and 
B coefficients denoted as acoeff and bcoeff, which will be discussed later in this section. The A and B 
coefficients along with a and fi factors are used in the calculations of K t and K p gains. 

For each Specld test, a maximum Specld temperature, a minimum Specld temperature, and each zone low 
power level per cycle (not per second) must be preplanned. The maximum Specld temperature should not be 
too high since it may damage a test specimen in an open test environment, and the minimum Specld 
temperature should not be too low since it will take a long time to finish the test. There are two phases in a 
Specld test: a heating phase with low power level and a cooling phase with no power. When any control TC 




reaches the maximum Specld temperature; the heating phase is finished, power is turned off, and natural 
cooling occurs. When all control TCs reach the minimum Specld temperature, the Specld test is complete. 
New acoeff and bcoeff coefficients are defined as follows in equations (32) and (33): 


acoeff = —a 

(32) 

1 

bcoeff - — 
b 

(33) 


Coefficient acoeff is calculated during the cooling phase. Consider equation (1 1): 

dc 

— — ac + bu + q (11) 


dc 

At t = 0, u = 0, — = 0 , and c = c 0 where c 0 is the initial ambient temperature and g - -ac 0 . Substitution of 
dt 

g = -ac 0 into equation (11) yields equation (34): 


dc 

dt 


— a(c — c 0 ) + bu 


(34) 


During the cooling phase there is no power; so u = 0 and as shown in equation (35): 

dc 


dt 


— a(c — c 0 ) - — a(c 0 — c) 


Replacing equation (32) into equation (35) can be expressed as equation (36): 

dc ( 1 


dc ( 1 \ 

acoeff - — ( 

dt Vc 0 — cj 


(35) 


(36) 


During the heating phase, u is a low constant power level. Using equations (32) and (33), u can be 
expressed in acoeff and bcoeff coefficients as equations (37) and (38): 


1 r dc 


u = 


dt 


~ a(c ~ c 0 ) 


bcoeff 


r dc 
dt 


+ acoeff(c — c 0 ) 


(37) 


bcoeff = -T- 

+ acoeff (c - c 0 ) 


(38) 


Equations (36) and (38) are used during Specld to calculate acoeff and bcoeff coefficients. For a real-time 
computer system, the calculation of power levels is carried out in each time n, and as mentioned earlier, F is 
the number of cycles per second. The power level at time n can be expressed as follows in equation (39): 
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1 r [dc 1 ) 

U n = p \bcoeff 1^— + acoeff (c - c 0 )Jj 


( 39 ) 


Equation (39) is the controller output equation derived from the heating process model and is only used for 
the first five seconds at the beginning of a thermal test for any QL heater zones. 

Figure 8 shows in detail how acoeff and bcoeff are calculated. During the cooling phase, C H2 is the starting 
temperature at the start time t H2 , and C L2 is the ending temperature at the end time t L2 of the selected 
temperature range, c c is the midpoint temperature of c H2 and c L2 as shown in equation (40): 


C H2 + C L2 



and the slope of skin temperature versus time is expressed in equation (41): 

dc _ c H2 — c L2 
dt t H2 — t L2 


(40) 


(41) 


Because of the delay, the temperature continues to rise for a short time after the heating phase is done, and 
the power is turned off. The first selected cooling temperature range should be chosen well below the 
temperature peak to ensure the equivalent heating temperature range is in the heating phase with u > 0.0 . The 
initial ambient temperature c 0 is always less than the midpoint temperature c c , where c 0 -c c < 0.0, and the 

temperature is always decreasing with time for the cooling phase, dc/dt < 0.0 ; therefore, acoeff is always 
positive. Equation (36) is used to calculate acoeff 


For the equivalent heating phase, the midpoint temperature c h must be selected as closely as possible to 
the cooling phase midpoint temperature C c within ±1.0 °F. Let C IA = c h — Ac be the starting temperature at 
the start time t Ll , and c in — c h + Ac is the ending temperature at the end time t in of the selected temperature 
range, and the slope of skin temperature versus time is shown in equation (42): 


dc _ c H1 — c L1 
dt t H1 — t L1 


(42) 


The initial ambient temperature c 0 is always less than the midpoint temperature c h , c h -c 0 > 0.0 , and the 
temperature is always increased with time during the heating phase, dc/dt > 0.0; therefore, bcoeff is always 
positive. Equation (38) is used to calculate bcoeff. 

The above calculation is for one control TC and for one selected temperature range. Many more selected 
temperature ranges should be used to calculate acoeff and bcoeff The final acoeff and bcoeff values are the 
average of all the calculated values. 
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The integral gain K j and the proportional gain K p can be determined from equations (15), (26), and (27). 
With the time interval At = 1 s, equation (15) becomes: B = b, and equation (27) becomes equation (43): 


K„ 


/? = /? 

B b 


(43) 


By replacing 1/ B = 1/ b = bcoeff into the above equation, K p can be expressed as equation (44): 


K p = fibcoeff (44) 

Equation (44) is the equation used in the calculation of K p . 

With the time interval At equal to one second, equation (26) can be written as equation (45): 

Ki = aK p (45) 

Equation (45) is the equation used in the calculation of if, . 

Power Control Cabinet 

In order to implement the TCS, an understanding of the operation of the PCC is necessary. The PCC is 
actually the main component of the heating process mentioned in the QL TCS design. For each control TC at 
time n, the ACP uses the delayed power level U n _ d , the skin temperatures C n _ d and C n , the thermal setpoints 
R n _ d and R n , K i , and K p to calculate power level U n which is a real number ranging from 0 to 60/ F power 
levels. Since the ACP calculates the power level 5 times per second (F= 5), U n value ranges from 0.0 to 12.0 
power levels per cycle. The ACP translates U n into a format suitable for interfacing with the output module 
or D/A converter. After receiving the signal from the ACP, the output module converts it into the current 
signal / between 4.0 mA and 20.0 mA and sends it to the PCC which fires the QL heater installed above the 
control. TC attached on a test article or an individual test coupon. Then, the temperature feedback from the 
control TC is input into the ACP for closed-loop control. 

The PCC is a back-to-back silicon controlled rectifier (SCR) zero crossing controller that operates on 
60 Hz or 60 power levels per second (ref. 5). For the circuit configuration in figure 9, the SCRs must be 
connected in parallel, but opposite directions. The forward SCR conducts during the positive half of the cycle, 
and the reverse SCR conducts during the negative half of the cycle to allow both halves of the AC waveform. 
In AC, the zero-crossing is the instantaneous point with no voltage or current where the sine wave goes from 
negative to positive or from positive to negative. The zero crossing detector senses the zero crossing of the 
AC voltage/current input, and the SCR is turned on. The benefits of zero crossing are very low harmonics and 
radio frequency interference. The PCC has zero crossover firing mode that distributes even on/off cycles to 
provide the desired percentage output based on the current signals in mA sending to the SCRs. The QL heater 
is fired when the zero-crossing occurs where the sine wave goes from negative to positive as shown in figure 
10a. Table 1 and figures 10a to 10c demonstrate the relationships between the current signals /to the SCRs, 
percentages of full power, number of cycles on/off, power level per second, and power level per cycle. 
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Table 1. Relationships between current signals / to SCRs, percentage of full power, number of cycles on/off, 
and power levels. 


Signals to SCRs, 
mA 

Percentage of 
full power 

Number of cycles 
on/off 

Power levels 
per second 

Power levels 
per cycle 

20.0 

Full on 

all cycles on 

60 

60 IF 

16.0 

75% 

3 cycles on, 
1 cycle off 

45 

45 IF 

12.0 

50% 

1 cycle on, 
1 cycle off 

30 

30 IF 

8.0 

25% 

1 cycle on, 
3 cycles off 

15 

15 IF 

4.0 

Full off 

All cycles off 

0 

0 


QUARTZ LAMP THERMAL CONTROL SYSTEM IMPLEMENTATION 

Equations (9) and (39) by themselves cannot build the TCS; other important information is required to 
build the TCS. The required information such as test type, heater type, zone type, and thermal profiles was 
preplanned, created, compiled, and loaded in binary format in the ACP before any thermal test and is called a 
thermal map. A thermal map contains a series of important information that instructs the TCS to perform its 
intended tasks. The implementation of the QL thermal control algorithm written in ANSI C was based on the 
derived thermal control equations and thermal maps. 

To implement thermal control, besides coding, a small four-zone setup was built with 4 input modules, 
4 Type-K control TCs, 4 individual test coupons, 4 output modules, and 4 infrared heated lamps. The test 
setup was covered with an insulation blanket to minimize the effect of natural and forced convection. The test 
setup was limited to a maximum temperature of 350 °F for safety. When the TCS algorithm proved to be 
working with the small setup, a larger setup similar to the one shown in figure 1 1 with 18 QL zones and 18 
different test coupons separated by insulated ceramic walls. Each zone has a QL heater installed on the lid 
and a Type-K control TC attached on top of a test coupon. Figure 12 shows a typical control TC attached on 
top of a test coupon. Figure 13 shows a typical single quartz lamp bulb. Figure 14 shows a QL carbon/carbon 
(C/C) elevon large-scale heating test in the FLL using the PI TCS. 
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Thermal Control Equations 

The key of the thermal control system is the derived thermal control equations (9) and (39) from which 
the ACP calculates power level every 200 ms. Equation (9) is the closed-loop control equation and the main 
equation to calculate power level. 

Equation derived from the PI controller 

Power level at time n is calculated using the closed-loop thermal control equation (9): 



(9) 


The calculated power level ranging from 0.0 to 12.0 per cycle must be converted to a suitable integer number 
to the assigned output module (D/A). Upon receiving the input signal, the output module sends the 
corresponding current signal to the assigned PCC channel which will fire its connected QL heater. The zone 
control TC sends the feedback temperature to the ACP via the assigned input module (A/D). The ACP uses 
the feedback temperatures along with the profile setpoints and the power level at time n-cl to calculate the 
power level at the current time n. 

Equation derived from the heating process 

Power level at time n is calculated using equation (39): 


U n = bcoeff 


dc 

— + acoeff(c - c 0 ) 


(39) 


At the beginning of a thermal test, equation (39) is only used for 5 s to calculate power level at time n. When 
enough data is saved, equation (39) is no longer used. 


Thermal Maps 

Three thermal maps: a profile control map (PCM), a closed-loop map (CLM), and a power output map 
were created to implement the TCS. Each map contains important information that is used in the thermal 
control algorithm and is described in the following sections. 

Profde Control Map 

The PCM consists of preplanned thermal profiles in segment format for all control TCs. Each segment 
contains start time, end time, start temperature, end temperature, and other important information. The ACP 
uses the data programmed in each segment to calculate thermal setpoints and to make decisions in executing 
its tasks. The programmed capabilities such as limit checking, average power level, maximum power level, 
and feedback temperature sensor switching are preplanned in each segment. The graph in figure 15 displays 
a thermal profile defined in Appendix E that has 8 segments with the highest temperature at 1480 °F and the 
thermal profile time of 1200 s. 

The PCM has the number of control TCs, test type (Specld, replay power level, closed-loop, or open-loop 
power level), heater type (QL, or GH, or QL GH), PCC used in the test, number of control delay cycles, 
Specld maximum temperature, Specld minimum temperature, maximum ramp up rate, and more. 


13 



Closed-Loop Map 

The CLM holds the information that links which PCC channel and which control TC belong to which 
zone, which zone is master, and which zone is slave. A master zone always has its assigned control TC, and 
power level is calculated for this zone using the thermal control equation (9). A master zone can have none, 
one, or multiple slave zones. A slave zone does not have any control TC and is assigned to a master zone. A 
slave zone is fired with the same power level as its master zone. The CLM also has each zone heater type (QL 
or GH), Specld power level, temperature limits including highest, lowest, high deviation, and low deviation. 
A deviation equals to a thermal setpoint at time n-d R n _ d minus its corresponding skin temperature at current 
time C n (deviation = R n .d - C„). The TCS is underheated when the deviation is positive and overheated when 
the deviation is negative. There are three zone configurations: one PCC channel for one control TC (most 
common), one PCC channel for multiple control TCs (feedback temperature sensor switching), and one PCC 
channel with no control TC (slave zone). Since a zone can have one or multiple control TCs, and also since 
the slave zone doesn’t have a control TC, the number of PCC channels and the number of control TCs may 
not be the same. 

Power Output Map 

The PCM contains the number of zones and the physical addresses of all output modules used in the test. 
The zone defined in the CLM is linked to this map. 

THERMAL CONTROL SYSTEM CAPABILITIES 

The TCS capabilities are categorized into three different types: interactive, programmed, and pre-test. The 
interactive type must be triggered by an operator’s action. The programmed type must be preplanned and 
programmed in the thermal maps. The pre-test type must be done before actual tests. 

Interactive Capabilities 

The TCS provides the capabilities to interactively control a thermal test (ref. 3) by the operator such as 
ramp up, start, hold, resume, and stop. The operator can initiate a thermal test by selecting the ramp up or start 
function. The operator can hold a thermal test by selecting the hold function and can resume a thermal test by 
selecting the resume function. The operator can stop a thermal test at any time by selecting the stop function. 

Ramp Up 

During the ramp up, the TCS brings all zones to their first thermal setpoints and maintains these 
temperatures while waiting for thermal start to be initiated. For each active control TC, the ramp up 
temperature difference between its first thermal setpoint and its initial skin temperature is computed. The 
maximum ramp up temperature difference is determined by comparing all the ramp up temperature 
differences. The ramp up time is calculated through dividing the maximum ramp up temperature difference 
by the maximum ramp up rate defined in the PCM. From the ramp up time, the ramp up rate for each control 
TC is calculated via dividing the ramp up temperature difference by the ramp up time. In order for the first 
thermal setpoints to start at the Thermal Profile Elapsed Time (THRM PET) equal to zero, the ramp up time 
must be a negative number as displayed in the real-time temperature graphs. When the ramp up finishes, the 
TCS will inform the operator and maintain the first thermal setpoint for each zone at THRM PET equal to 
zero. 

Start 

If the ramp up had been selected and the TCS finished the ramp up, the TCS would start segment 1 
immediately when the start function is selected. If only the start is selected, the TCS will go through the ramp 
up period described above. When the ramp up time is reached, the TCS starts segment 1 at THRM PET equal 
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to 0. For each control TC and each segment, the TCS uses the information in the PCM to compute thermal 
setpoints and to execute its tasks. When one segment ends, the next segment starts until the last segment is 
completed. During a thermal test, the TCS calculates the power level for all control TC every 200 ms until 
THRM PET is equal to the preplanned thermal profile time programmed in the PCM. When a thermal test 
finishes successfully, the TCS zeroes out all power level buffers, stops the thermal test, shut down all PCCs 
used in the test, and notifies the operator. 

Hold 

After the ramp up or the start is selected, the hold function becomes available. When the hold is selected, 
the TCS holds each control TC profile temperature at the thermal hold time, and THRM PET is not increased 
until the resume function is selected. 

Resume 

When the resume function is selected, the TCS resumes from where the hold function has been selected. 
THRM PET increases, and the profile temperatures continue from the thermal hold time. 

Stop 

After the start or the ramp up function is selected, the stop function can be selected at any time during a 
thermal test. When the stop function is selected, the TCS zeroes out all power level buffers, stops the thermal 
test, shuts down all PCCs used in the test, and notifies the operator. There will be no firings and no options to 
go back to the thermal test from this point. 


Programmed Capabilities 

There are additional TCS capabilities that are programmed in the thermal maps: limit checking, average 
power level, maximum power level, feedback temperature sensor switching, and replay power level. 

Limit Checking 

The limit checking capability offers the option to shut down a thermal test in case any control TC 
temperature exceeds its predefined limits programmed in the CLM. There are four limits: highest temperature 
limit, lowest temperature limit, high deviation limit, and low deviation limit. The limits can be defined with 
any appropriate values based on each test. The highest temperature limit and the high deviation limit are for 
the overheated cases. The lowest temperature limit and the low deviation limit are for the underheated cases. 
In each segment, the information programmed in the PCM provides the TCS with which types of limit 
checking that must be performed. 

Average Power Level 

The average power level capability provides the TCS a choice to either use the calculated power level or 
average it with previous power levels. This capability prevents sudden changes in power levels. In the PCM, 
there is an integer representing the number of previous cycles used in the average power level. If this integer 
equals to 0, the calculated power level is used. If this integer is a non-zero, an average power level is computed 
using 50 percent of the calculated power level and 50 percent of previous power levels specified by the integer. 

Maximum Power Level 

The maximum power level capability allows putting a limit on the power level in each segment. If the 
calculated power level is greater than the programmed maximum power level in the PCM, it is reduced to the 
maximum power level. This capability helps to prevent overheating. 
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Feedback Temperature Sensors Switching 

The feedback temperature sensor switching capability provides switching from one feedback temperature 
sensor to another feedback temperature sensor in one zone. For example, with the thermal profile in Appendix 
E, zone 1 has two control TCs, CTC01 and CTC02. CTC01 is the feedback temperature sensor from segment 
1 to segment 4. CTC02 is the feedback temperature sensor from segment 4 to segment 8. The TCS checks 
every segment to know when to switch the feedback temperature sensor. 

Replay Power Level 

The replay power level capability allows running a thermal test with power levels recorded in a previous 
thermal test, and either continuing to closed-loop thermal control or finishing the test. The selected power 
level recorded file must be loaded into the ACP before starting the test. In the PCM, the switch time and the 
transition time must be programmed. When the thermal test starts, the power level recorded file is opened, 
and the power level data is read until non-zero power levels are found. Every cycle from here, the TCS 
compares THRM PET with the switch time. If THRM PET is less than the switch time, the TCS reads power 
levels from the recorded file and fires the heaters. After the switch time is reached, the TCS checks for the 
transition time. If the transition time equals to zero, the thermal test is done, and the TCS zeroes out all power 
level buffers, stops the thermal test, shuts down all PCCs used in the test, and notifies the operator. If the 
transition time is greater than zero, the TCS performs a transition from replay power level to closed-loop 
control and brings up all zones to their first profile temperatures during the transition time. After the transition 
period, the TCS automatically starts segment 1 of the thermal profiles and continues until the last segment. 

Pre-test Capabilities 

After a thermal test setup is completed, some preliminary thermal tests must be conducted before running 
actual tests. Preliminary tests include lamp check, zone check, and Specld. These tests help to avoid 
unnecessary failures during actual tests and contribute in saving time, resources, money, and maintaining 
on-time schedules. 

Lamp Check 

The lamp check capability allows checking each individual output module (D/A) that is connected to its 
PCC channel and QL heater. Selecting an output module and selecting a current signal from 4 m A to 20 m A 
(normally 5 mA is used), an operator can fire the corresponding QL heater. If the corresponding QL heater 
doesn’t fire as instructed, there must be a problem with either the output module or the associated PCC channel 
or the QL heater and the problem must be corrected. 

Zone Check 

The zone check capability is used to verify the setup of each zone with its output module, PCC channel, 
QL heater, and control TC. Each zone is commanded with a very low power level in the range of 3 percent to 
10 percent of full power; normally 3 percent is used. Test personnel with proper eye protection performs a 
visual check and verifies that the correct PCC channel and QL heater are fired. The zone control TC is checked 
and verified for rising temperature. This capability allows zones to be checked in any combination from one 
at a time to all zones at once. 

Specimen Identification 

After the lamp check and the zone check, a Specld test is run to determine acoeff and bcoeff coefficients. 
The Specld was explained in detail in the integral gain coefficient, K ; , and proportional gain coefficient, 
K , section. A Specld maximum temperature, a Specld minimum temperature, and a very low Specld power 
level for each zone are programmed in the thermal maps. When the Specld test is complete, acoeff and bcoeff 
coefficients are calculated, recorded, and ready for use. These coefficients are used in the calculation of K t 
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and K . If a test article or coupons were used in previous tests, the Specld does not need to be run since the 
acoeff and bcoeff coefficients were already calculated. 


GRAPHITE HEATER CLOSED-LOOP THERMAL CONTROL 

In 2005, the FLL purchased graphite heater PCCs with the objective to support heating tests with 
temperatures up to 3000 °F. Each graphite heater PCC has a built-in TCS and can be used by itself to run only 
GH. On the other hand, the TCS can be used for only QL heaters. Therefore, to have both QL and GH 
synchronized in the same test, GH closed-loop thermal control was implemented in the TCS. 

In 2006, the small test chamber was setup with two GH zones. The same thermal profile defined in 
Appendix E and shown in figure 15 was applied to both zones. Figure 16 displays the GH elements used in 
thermal tests. At the beginning, the derived QL thermal control equations were implemented for GH the same 
way as for QL. However, the QL thermal control equations did not work for GH. The real-time deviation and 
power level recorded data were studied and analyzed. The real-time graph in figure 1 7 shows a long delay 
and the cyclic pattern of overheating and underheating. At THRM PET equals to 200 s, the profile temperature 
is 530 °F; the real-time graph in figure 17 shows the skin temperature is around 275 °F, 255 °F below the 
desired temperature. The test data confirm s the sluggishness of the GH at the beginning of a test. It is found 
that when the system was underheated, the TCS produced high power levels in an effort to try catching up 
with the profile temperature; and when the profile temperature was caught up, it was soon overheated. 
Similarly, when the system was overheated, the TCS produced low power levels in an effort to try descending 
down to the profile temperature; and when the profile temperature was achieved, it was soon underheated. 
Since the GH elements have very high thermal inertia as compared to QL, these problems were expected. 

The solution for the GH beginning sluggishness is to fire with high power levels for the first five seconds 
depending on the ramp up temperature difference between the first thermal setpoint and the ambient 
temperature. For the underheating problem, the current deviation must be in the range of 1 5 °F to 0 °F, and 
the deviations of the last three cycles must all be higher than the current deviation; the control logic decreases 
the power levels depending on the values of the current deviation. For the overheating problem, the current 
deviation must be in the range of -15 °F to 0 °F, and the deviations of the last three cycles must all be smaller 
than the current deviation; the control logic increases the power levels depending on the values of the current 
deviation. The extra software control logic was developed to handle the GH high thermal inertia. The real-time 
graph in figure 1 8 demonstrates that the GH closed-loop thermal control operates within the required range. 
At THRM PET equals to 200 s, the profile temperature is 530 °F; the real-time graph in figure 18 shows the 
skin temperature is around 530 °F, very close to the desired temperature. 

The test setup was changed to one GH zone and one QL zone with the same previous thermal profile 
defined in Appendix E. The test results proved that the TCS new capability of combining QL and GH in a 
same test was achieved as shown in figure 19. At profile elapsed time equals to 200 s, the profile temperature 
is 530 °F; the real-time graph in figure 19 shows the skin temperature is around 530 °F, very close to the 
desired temperature. This combination test proves that the TCS can run with both QL and GH heaters in a 
same test. 


GRAPHITE HEATER OPEN-LOOP POWER LEVEL 

Later, there was a need to have the GH open-loop power level capability to heat a test article with 
predefined power levels. With the new capability, the TCS can fire the GH with power levels from a time 
history power level profile predefined by the user in the PCM and with no ramp up period. The hold, resume, 
and maximum power level capabilities are required to be the same as the closed-loop thermal control. The 
control TC for each zone must also be installed for monitoring the zone temperature. 
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In the PCM, power level profiles replace temperature profiles. In each segment, there are start time, end 
time, slope in power levels/sec, intercept in power level, and other important information. For this test type, 
Specld doesn’t need to be run. Only the highest temperature limit check is available and must be enabled to 
ensure that the open-loop power level test doesn’t run away from the predefined highest temperature limit, 
and bum the test article and/or the test setup. 

During a test if the hold function is selected, each zone continues to fire with the same power level at the 
thermal hold time, and THRM PET is not increased. When the resume function is selected, THRM PET is 
increased, and the power level profiles continue from the thermal hold time. 

A small test chamber with one zone of GH was used to develop this new capability. Verification and 
validation for the open-loop power level capability are accomplished by verifying the power levels in the 
recorded file with the predefined time history power level profile. 

The GH open-loop power level capability was utilized to heat a ceramic matrix composite (CMC) panel 
with preplanned power levels programmed in the PCM. The hold function was selected during a test to fire 
the GH at a certain power level for a period of time before returning back to the power level profile. There 
was also a test where the CMC panel was heated using the maximum power level of 12. All of these tests 
were conducted to study how fast the temperatures rise for different power levels. The highest limit check is 
always used in these tests to ensure that the test article and/or the test setup will not be overheated above the 
permitted temperature. When the temperature reached its highest temperature limit, the TCS would stop the 
test immediately to prevent any damages from occurring. 

LARGE-SCALE HEATING TESTS 

For almost a decade, this in-house PI TCS has fulfilled its responsibility to support large-scale heating 
tests in the FLL. All large-scale tests were accomplished in a large test chamber with high number of zones 
and large test articles. The TCS successfully met the test objectives defined for each test. The test data were 
used by researchers and engineers to evaluate the performance of the hot structures in the tests. 

The C/C control surface tests were the first QL large-scale heating tests that used the PI TCS (fig. 14). 
The carbon-silicon carbide (C/SiC) hot structure subcomponent tests were soon followed successfully. The 
above tests were the first combined thermal/mechanical tests performed on CMC hot structures. A series of 
QL thermal and thermaPmechanical tests were completed on other hot structures which included both C/SiC 
and C/C subcomponents. The combined QL GH thermal tests were executed on a test article where some GH 
zones achieved temperatures up to 3000 °F; the highest temperature ever accomplished by the PI TCS. 

THERMAL CONTROL SYSTEM DISCUSSION 

In this section, real-time graphs and test data of a typical CMC hot structure test are used to highlight the 
performance of the TCS. This thermal test has 18 QL zones, and the thermal profile time programmed in the 
PCM is 1900 s. Zones 2, 10, 12, and 13 are selected for the discussion. 

THRM PET in figure 20 displays 1706 s, 194 s short of the thermal profile time programmed in the PCM. 
The termination was caused by the high deviation limit of CTC34 (Zone 10) at THRM PET equal to 1706 s. 
The TCS stopped the thermal test, zeroed out all power level buffers, shutdown all PCCs used in the test, and 
informed the operator of the cause of the shutdown. The thermal test had stopped, but the test data recording 
went on until the test article cooled down to a safe temperature range; therefore, the real-time graphs 
continued. After the test stopped, the deviation and the power level recorded files were downloaded and 
analyzed. From the deviation data as shown in figure 21a, the shutdown of high deviation limit was verified 
and validated. At THRM PET equal to 1706 s and test time equal to 6622 s (fig. 20), CTC34 absolute value 
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of temperature deviation was 300.90 °F, just over the predefined high-deviation limit of 300.0 °F in the CLM. 
The TCS shutdown the thermal test instantly at the first time the high- deviation limit occurs. From the power 
level data shown in figure 21b, the power level went to zero at THRM PET equal to 1 600 s and test time equal 
to 65 1 6 s (fig. 20) which was 1 06 s before shutdown. Therefore, even without power for the last 1 06 s before 
shutdown, the CTC34 remained at much higher temperatures than the preplanned profile temperatures and 
caused the high-deviation limit shutdown. Overheating during the cool down phase with negative slopes 
(-5.81 °F/s for CTC33 and -5.31 °F/s for CTC34) happens often in test setups not incoiporate with an active 
cooling system. Actually many other control TCs would exceed high-deviation limits if CTC34 did not. 

After the thermal test is terminated, figure 22 shows that the thermocouple junction of CTC39 (Zone 2) 
opened, and its temperature jumped over 3200 °F at test time equal to 6675 s. This real-time display provided 
the test team visual information that CTC39 must be fixed before the next test. 

CTC15 (Zone 12) and CTC16 (Zone 13) had the same thermal profile. Zone 12 and Zone 13 real-time 
graphs in figure 23 were overlapped, and it is impossible to distinguish them during the heating phase. Only 
during the cooling phase, the two graphs started to separate at THRM PET equal to 1630 s and test time equal 
to 6546 s. The graphs of the preplanned thermal profiles for CTCl 4, CTCl 5, and CTCl 6 are shown in figure 
24. In order to compare the real-time graphs and the preplanned thermal profiles, the ramp up time and the 
first segment start time are mentioned. From the test data, ramp up started at THRM PET equal to -43 s and 
test time equal to 4873 s, and segment 1 started at THRM PET equal to 0 s and test time equal to 4916 s 
(fig. 23). The graphs of real-time deviation and power level for Zone 12 are in figure 25, and the graphs of 
real-time deviation and power level for Zone 13 are in figure 26. Even Zone 12 and Zone 13 real-time 
temperature graphs were inseparable, the deviations and power levels were different at all times. For Zone 12, 
the CTCl 5 deviations as shown in figure 25a were still in good operating range, and Zone 12 power levels as 
shown in figure 25b were very small but non-zeroes at THRM PET equal to 1630 s and test time equal to 
6546 s (fig. 23). For Zone 13, the CTC16 deviations as shown in figure 26a started to be negative and 
continued increasing more negative, and Zone 13 power levels as shown in figure 26b had already been zeroes 
at THRM PET equal to 1630 s. The above discussion explains why Zone 12 and Zone 13 real-time graphs 
started to separate from each other during the cooling phase with no power present. 

The deviation and power level graphs demonstrate that at the beginning, the QL thermal control is very 
oscillatory at low temperatures and low power levels with more off cycles. When a quartz lamp is fired, its 
control TC temperature increases; and when the quartz lamp is off, its control TC temperature decreases. The 
control becomes much better at higher temperatures and higher power levels with more on cycles. Since this 
TCS does not have an active cooling system, one of the zones with negative slopes during the cooling phase 
normally triggers the shutdown due to exceeding high deviation limit. Even with no power, all zones retain 
so much heat causing the control TC temperatures to be too high compared to the corresponding profile 
temperatures. 

For all thermal tests, the construction of a zone, and the positions of control TC and heater are very 
important with respect to the performance of the TCS. If a zone is insulated well, and the control TC is 
positioned near the heater; the control is very good, and the errors can achieve ±0.5 °F. If a zone is insulated 
poorly, and/or the control TC is positioned not well; the control can become oscillatory, and the errors can 
reach ±20.0 °F. 


CONCLUDING REMARKS 

The PI thermal control system described in this paper is a very reliable, safe, and robust system that was 
built in-house to support thermal testing in the FLL. The derived QL thermal control equations and algorithm 
are implemented with ANSI C in the ACP, a component of the existing DACS. The ACP calculates power 


19 



levels every 200 ms and sends current signals to the PCCs to fire QL and/or GH heaters to heat up a test article 
following preplanned thermal profiles. 

From 2003 until 20 1 0, the PI thermal control system has performed many thermal tests as described in 
the Large-Scale Heating Tests section. There are significant advantages in having an in-house thermal control 
system. When new requirements come up, new software and new test setups can be developed, tested, verified, 
and validated quickly. The GH closed-loop thermal control and GH open-loop power level are added later to 
satisfy the demand of higher temperature tests (up to 3000 °F). Also, when a thermal test fails, the deviation 
and the power level recorded files can be used to analyze the causes of the problems, and the software fixes 
can be fast, helping the projects to stay on schedule. 

Theoretically, equation (9) is the general PI control equation that can be used to build any control system 
with R as the reference input, C as the controlled variable, E = R-C as the difference between R and 
C , and U as the controller output. An appropriate continuous-time mathematical model for the controlled 
process similar to the heating process must be defined between the controlled variable C and the controller 
output U . Then, the continuous-time controlled process equation must be transformed to the discrete-time 
controlled process equation. The U in the controlled process equation is replaced with the U in equation (9) 
to develop the relationship between C and R from which the transfer function C(z)/R(z) in the z-domain is 
derived. For the system stability, the poles of the transfer function should be assigned on the real axis in the 
right-hand side inside the unit circle (0.0 < zi, ... Zi < 1.0) to simplify the mathematical process. From the 
derivation of the transfer function and the controlled process equation, a stability region of the control system 
can be determined. The values of K i and K p should be tuned and selected carefully so that the control system 

is stable and robust. Besides the derived control equations, the control algorithm should always have 
capabilities similar to the ones described in the Thermal Control System Capabilities section such as reference 
input limit checking, and average and maximum controller output value to make the control system more safe 
and versatile. 
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FIGURES 
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Figure 1. A quartz lamp thermal control system block diagram. 
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Figure 2. Proportional integral controller block diagram. 
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Figure 3. The integral of e = r — c in A t = t n — t n _ t . 


Controlled output Controlled variable 



150138 


Figure 4. Fleating process model block diagram. 
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Figure 5. Closed-loop PI thermal control system block diagram. 



Figure 6. Selected control system stability region in the z-domain. 
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Figure 7. The PI thermal control system stability region. 



Figure 8. Determination acoeff and bcoeff coefficients during Specld. 
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Figure 9. The power controller back-to-back silicon controlled rectifier configuration. 



Figure 10a. 100 percent of full power firing, all cycles on. 



Figure 10b. 75 percent of full power firing, 3 cycles on and 1 cycle off. 

Figure 10. Relationship between percentages of full power and number of cycles on/off. 
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Figure 10c. 50 percent of full power, 1 cycle on and 1 cycle off. 
Figure 10. Concluded. 
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Figure 1 1. A multi-zone quartz lamp heater test setup. 
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Figure 12. A typical control thermocouple attached on top of a test coupon. 
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Figure 13. A typical single quartz lamp bulb. 
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Temperature, 



150150 


Figure 14. A QL C/C elevon large-scale heating test in the FLL. 



Figure 15. Graph of the thermal profile described in Appendix E. 
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150152 


Figure 16. Graphite heater elements. 
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Figure 17. Real-time temperature graphs of two GFI zones. 
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Figure 18. Real-time temperature graphs of two GH zones with control logic software. 
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Figure 19. Real-time temperature graphs of one GFI zone and one QL zone. 
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Figure 20. Real-time temperature graphs of Zone 9 and Zone 10. 
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Figure 21a. Real-time graph of temperature deviation. 

Figure 21. Real-time graphs of temperature deviation and power level for CTC34, Zone 10. 
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Figure 21b. Real-time graph of power level. 


Figure 21. Concluded. 



Figure 22. Real-time temperature graphs of Zone 1, Zone 2, Zone 3, and Zone 4. 
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Figure 23. Real-time temperature graphs of Zone 11, Zone 12, and Zone 13. 



Figure 24. Preplanned thermal profiles of Zone 11, Zone 12, and Zone 13. 
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Figure 25 a. Real-time graph of temperature deviation. 
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Figure 25b. Real-time graph of power level. 

Figure 25. Real-time graphs of temperature deviation and power level for CTC15, Zone 12. 
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Figure 26a. Real-time graph of temperature deviation. 
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Figure 26b. Real-time graph of power level. 

Figure 26. Real-time graphs of temperature deviation and power level for CTC16, Zone 13. 
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APPENDIX A 

Derivation of PI Controller Output 

In this appendix, the discrete PI control equation (6) is written as equation (Al): 

U„ +K,y(£„ +E„_,) + K p (E„ 

At time n = 0, equation (Al) becomes equation (A2): 

t/ 0 (U 0 = 0) 

At time n = 1, equation (Al) becomes equation (A3): 

[/,=[/„ + K, -j(E, +£„) + K p (E , - £„) 

Equation (A3) can be re-written as equation (A4): 

U, =£/„ +X,y[2(£, + £„)-£, -£„]+£,(£, -£ 0 ) 

At time « = 2, equation (Al) becomes equation (A5): 

u 2 = u, + K A (£ 2 + £, ) + Ar p ( E 2 - £ , ) 

Replace C j in equation (A4) into equation (A5), and equation (A6) is achieved: 


U 2 =U 0 +K,^[2(E I+ E 0 )-E t -E 0 ] + K p (E,-E 0 ) + 

£,y(£; 2 +E l ) + K p (E 2 -E l ) 

Expand and combine the like terms in equation (A6), and equations (A7) and (A8) are derived: 

U 2 =U 0 + K^[2E X +2E 0 -E l -E 0 + E 2 +E l ] + 

K p (E l -E 0 + E 2 -E l ) 

U 2 =U 0 +K i ^-[2(E 2 +E X + E 0 )-E 2 -E 0 ] + K p (E 2 -E 0 ) 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 
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At time n = 3, equation (Al) becomes equation ( A9): 


u, =u 2 + K i -(E 3 + E 2 ) + K p (E 3 -E 2 ) 


Replace U 2 in equation (A8) into equation (A9), and equation (A 10) is achieved: 


U, =U„ +K,-[2(E 2 +E, +E,)-E 2 -E„]+K p (E 2 -EJ + 
kE(E 2 +E 2 ) + K p (E 2 -E 2 ) 

Expand and combine the like terms in equation (A10), and equations (Al 1) and (A12) are derived: 

u, = u„ + K, ^-[2(E 2 + E l+ E„)-E 2 -E„ + E 2 + E 2 ]+ 

K (E 2 — E 0 + E 3 — E 2 ) 


(A10) 


(All) 


U 3 =U 0 + K—[2(E 3 +E 2 +E x +E 0 )-E 3 -E 0 ] + 
K p (E 3 -E 0 ) 


(A 12) 


From the above expansion of equation (Al) for n = 0, 1, 2, 3, a general equation for the controller output 
power level U n at time n is derived and shown in equation (A13): 


U n =U 0 + K i — 2^Ej -E n -E 0 +K p (E n -E 0 ) 


(A 13) 


Moving the factor — into the bracket [ ] and combining E n and E 0 term, equation (A14) is obtained: 

U„ = U 0 +K : At ^£.- ( £ .+ £ o) +K p < £,-£,) 

j = 0 2 


(A 14) 


Equation (A14) is the form of the power level at time n for the time interval At. 
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APPENDIX B 

Relationship between Controlled Variable and Reference Input 

In this appendix, the discrete heating process equation (21) is written as equation (Bl). 


c„ + ,= 4c„-c 0 Mt/„+c 0 


(Bl) 


A t in the heating process was assumed to be the same as A t in the control model. 
The discrete PI control equation (A13) was developed in Appendix A. 


At 

U n =U 0 +K i — 


2 H E j 

V AO J 


■E„ -E 


+ K(E n -E,) 


(A 13) 


Replace equation (A13) above into equation (Bl), and equation (B2) is achieved: 

V „ V 


C H+X =A(C n -C 0 ) + BUo+BK,* 


2 Z E J 


-E n -E o 


V AO J 
BK p (E n -E 0 ) + C 0 


+ 


(B2) 


Replace U 0 = 0 and E n =R n —C n in equation (B2), and equation (B3) is derived: 


C n +\ ~ AC n -AC 0 + 

r „ 


+ BK,At 




At 

~2 


Vi=o j = o J 

At 

2 


- BK — R + BK —C 

1 2 2 


-BK,—R 0 +BK i —C 0 


+ BK p R n -BK p C n -BK p R 0 +BK p C 0 

+ Cn 


Combine the similar terms of C„+;, C„, and R„ in equation (B3), and equation (B4) is obtained: 


(B3) 


r - r 

'-'« + ! — '-'M 


A + BK, — - BK n 

' 2 p 


f 


+ R, 


BK -BK, — 


V 


J 


f 


+ BKAt 




f 


\j = 0 AO J 

At 


At 


C 0 - AC 0 - BK, -R 0 + BK, — C 0 - BK p R 0 + BK p C 0 


(B4) 
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(B5) 


(B6) 


Equation (B6) demonstrates the relationship between the skin temperature C and the thermal setpoint R in the 
time domain. 
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APPENDIX C 

Derivation of Thermal Control Transfer Function 


In order to derive the transfer function of the thermal control system, the z-transform of discrete equation (22) 
is taken. In this appendix, equation (22) is written as equation (Cl): 


C n+1 = C„ 


A-BK n +BK 

p 1 2 


r 


+ R, 


BK-BK 

p 1 2 


r _n n ^ 

BK,At 2k “SC 

v 7= 0 7=0 J 


(Cl) 


+ 1 (\ 


The z-transforms for the two terms and X C' in equation (Cl) are complicated; therefore, extra 


7=0 


7=0 


intermediate steps are needed to make it easy to understand. The term X R ,■ is expanded as follows in equation 


7=0 


(C2): 


-^0 +^1 + ••• + R n 

7=0 


(C2) 


n 

From equation (C2), the term X R , can be represented graphically in figure Cl below. 

7=0 



Figure Cl. Graphical representation of X^, • 


7=0 


Let t/(/) be the unit step function occurring at / = jAt as shown in equation (C3): 


u(t - jAt) = 1 at t = jAt 


(j = 1, 2, ..., h) 


(C3) 
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'ZR: can be represented as a discrete convolution of r(t)and u(t) as demonstrated in equations (C4) and 

7=0 

(C5): 


n 

y Rj = r(o)u(t - nAt)+ r(\)u(t - (n - l)Af) + + r(n)u( 0) 


(C4) 


7=0 


n n 


IX =Yj r ^ u ^ n ~j) 

7=0 7=0 


(C5) 


n 

Similarly, £C - can be represented as a discrete convolution of c(/)and 7/ (/ ) as expressed in equation (C6): 


7=0 


n n 


IX = H c U) u ( n ~j) 

7=0 7=0 


(C6) 


By using the discrete convolution technique (ref. 2), the z- transforms of X R , and £C are as follows in 

7=o 7=0 


equations (C7) and (C8): 


z('£R j ) = R(z ) | 
j-o 


f z ^ 


yz-ly 


1 — z 1 


(C7) 


Z(Z C 7') = C (Z)| 

7=0 


f „ \ 


vz-ly 


C(z) 

1 — z 1 


(C8) 


The z-transform of C„ +1 can be written as follows in equation (C9): 


4c„ +1 ) = z(c(z) - c 0 ) 


(C9) 


Take z-transform of equation (Cl), and equation (CIO) is achieved: 


z(C(z)-C 0 ) = C(z) 


+BK: — 

p « 2 


i?(z) 


BK-BK — 

p * 2 


(CIO) 




V 


/?(*) C(z) 
1-z' 1 1 — z 1 


1 - z" 
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Replace z with — j- into equation (CIO), and equation (Cl 1) is obtained: 

z 


-h(C(z)-C 0 )=C(z) 


R(z) 


f At' 

A - BK + BK — 

V P 2 j 

f 

BK - BK 

V P 2 

A m „\ \ 


+ 


Jr 




5^0 r 


i?(z) C(z) 


Vl-z' 1 1 - z A j 


1 -z' 


(Cll) 


Multiply the left-hand side and the right-hand side terms with z' 1 or (l - z 4 ) or z A- z 4 ) so that both sides 
have the common denominator ofz -1 (l - z 4 ) to obtain equation (02): 


(l - z 4 )(C(z) - C 0 ) = z 1 (l - z- 1 )c(z)(j - BK p + BK, 
z- 1 (\-z- 1 )r(z) ( BK p -BK t 


At' 
2y 
At' 


( 02 ) 


J 


z- x BK t Al(R{z) - C(z)) + z"‘/ 0 

Expand the above equation and combine the alike terms, and equations (C13)-(C15) are obtained: 


C(z) 


1 — z 




At 

A-BK+BK 

p 1 2 


+ z l BK ; At 




R(z) 

z A [l 0 -C 0 ] + C 0 


At 

BK-BK,— 

p 1 2 


+ z BK-At 


+ 


(Cl 3) 


C(z) 


- z 1 (z 1 r 1 j A BK p + BK, — 


+ z BK.At 


J 


R(z) 
z- 1 (/ 0 -C 0 )+C 0 


( z 1 _ z‘ 2 )yBK p - BK t 


+ z BK ; At 


+ 


(Cl 4) 
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C(z) 


1 — z~ 1 + A-BK-BK- — 

{ p 2 




+ z 


A 

A-BK+BK 

p 1 2 




BK+BK,— 

p 1 2 


A A 
-2 
Z 


bk-bkA 


+ 


V 


(Cl 5) 


z- 1 (/ 0 -C 0 )+C 0 

Equation (Cl 5) is the z-transform of the skin temperature C(z) and the thermal setpoint R(z). Equation (Cl 5) 
will be used to derive the thermal control transfer function. 
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APPENDIX D 

Determination of Thermal Control System Stability 

To determine the TCS stability, the relationships between the integral sensitivity factor a, the proportional 
sensitivity factor p, the poles zi, and zi, in equation (Dl) and equation (D2) must be solved where 
a,j3 > 0, and 0 < z x , z 2 < 1 , therefore: 


0 < Zj,z 2 = 


2 - /? - —1 ± , f 2 - — T - 4 fl - — 

2 J VI ^ 2 J l H 2 


2-/?-— 1 - 4 fl -/?+— 1>0 
2 J l 2 J 


Stability Region 1 


Multiply equation (Dl) by 2 to obtain equation (D3): 


0 < (2-^-—) ±,1(2 -/?-—) — 4 f 1 — + — ) < 2 

2 t \ H 2 J l r 2 


Equation (D3) has two-side conditions, so consider the left-hand side in equation (D4): 


0 < r 2 — /? — —i + . {2- p ~ — 1 ~a{\~ p +— 
l 2 J vl H 2 J 1^2 


Move the square root term in equation (D4) to one side to achieve equation (D5): 


1 (2-/?-^T- 4 (l- ( 9 + S»)<2- ( 9-^ 
V v 2/1 2/ 2 


Square both sides as shown in equation (D6): 


2 _^Y 

2 J 2 J l P J 


Remove the |^2 - ft — — J term from both sides in equation (D6), and equations (D7) and (D8) are obtained: 

/ , D 

4 \- [i +— >0 (D7) 

v 2 
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Move all / 3 terms in equation (D8) to the left hand side as shown in equations (D9) and (DIO): 


4 1 


a 

2 


<1 


(D9) 


0<^ 


2 

< 

2 - a 


(DIO) 


To satisfy the conditions 0 < p in equation (D 1 0), the term 2 - a in the denominator must be positive as 
expressed in equation (Dll): 


0 < a < 2 (Dll) 

The boundary conditions of a in equation (Dll) will be applied to all stability regions. Equations (DIO) and 
(Dll) determine Stability Region 1 . 


Consider equation (D2): 


Stability Region 2 


( 2 -P- 


a/f\ 2 

2, 


2 


>0 


(D2) 


oB 

From equation (D8), 1 - ft + — > 0 ; therefore, the following condition shown in equations (D12) will always 
be true: 


abs\ 


a[l 




aft 


< abs\ 


2- p~— 
2 


(D12) 


aB 

if 2 -p~ ^ ^ 0 , then equation can be expressed m equation (D13^. 


■ abs. 


2 - [( - — 
2 


■41 1-/1 + — ) >{ 2 - p - 

' 2 ) \ H 2 


(D13) 


From the condition of equation (D13), both sides are negative and the left-hand side is less negative than the 
right-hand side. Therefore, equation (Dl) will have both poles z,, z 2 < 0 . Based on the selected stability 
condition that both poles must be greater than zero, the condition in equation (D 13) is not acceptable. 

(xB 

If 2 - p — — > 0 , then equation (D12) can be expressed in equation (D14): 



7 n a P\ a( -i n a P^\ (, n a P 

abs . 2-R- — -4 1-B + — < 2-B — — 

n 2 { ^ 2 J l H 2 


(D14) 


From the condition of equation (D14), both sides are positive and the left-hand side is less positive than the 
right-hand side. Therefore, equation (Dl) will have both poles z u z 2 > 0 which satisfy the left-hand side of 
the selected stability condition as shown in equation (D15): 


2 -/?- — > 0 
2 


(D15) 


Move all |3 terms in equation (D15) to one side, equations (D16) and (D17) are derived: 


2>/3 1 + 


(D 1 6) 


0 <^< 


2 + a 


(D 1 7) 


The boundary conditions of 0 < a < 2 in equation (Dll) produces the boundary conditions of ft as shown in 
equation (D18): 


0 < /? < 2 


(D18) 


Check if the condition 2 - f> — — > 0 in equation (D15) satisfies the right-hand side of equation (Dl): 


2- (j - — ) + ( 2 — /? — — Y — 4 ( 1 — /? + — ) < 2 

2 J vl 2 J l 2 J 


Move the 2 - // — — term in equation (Dl) to the right-hand as shown in equation (D19): 


+ ( 2 - /? - — ) - 4 f 1 - + — ) <2-f 2 -/7 - — 

VI 2 J 2 J 2 


Square both sides of equation (D19), and equation (D20) is derived: 


_ _ a/3\ ( a.p\ ( a/3') ( a/3\ 

2-y 6 — - -4 1 - 7 6+— < 4-4 2-t 8-— + 2-B — - 

2 l 2 J { ^ 2 { ^ 2 


Eliminate 2 - (3 term and divide by 4, and equation (D2 1 ) is obtained: 

V 2 1 


(D20) 
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(D21) 



<1- 


v 2 j 


l 2 J 


cxB 

From equation (D15), 2- [j - — > 0 ; therefore, the following condition in equations (D22) and (D23) will 


always be true: 


1- 


2-P-^' 
2 , 


<1 


f 

V 




a/3 


>-l 


(D22) 


(D23) 


The conditions in equations (D22) and (D23) satisfy the right-hand side condition of equation (Dl). Therefore, 
(xB 

the condition 2- ft- — > 0 satisfies both left-hand and right-hand sides of equation (Dl). Equations (D17) 
and ( D 18) determine Stability Region 2. 


Stability Region 3 

Expand equation (D2) as shown in equation (D24): 


2 p2 

,2 , « P 


4 + /3 z +—^ 4/3- 2a(} + ap - 4 + 4/? - 2a/? > 0 


(D24) 


Eliminate terms 4 and 4/3 and rearrange the similar terms in equation (D24), and equation (D25) is obtained: 


P 


f a 2 [i 


+ aft + ft -4a 


>0 


(D25) 


Since J3> 0, the term inside the bracket must be > 0 as shown in equation (D26): 


a 2 /? 


+ a/3 + /? - 4a > 0 


(D26) 


Move all /? terms in equation (D26) to one side, equation (D27) is derived: 


+4 a + 4)> 4 a 


(D27) 


Rewrite the above equation in a more compact form as shown in equation (D28): 
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— (a + 2) 2 >4 a 
4 v 7 


(D28) 


Move all a terms in equation (D28) to the right-hand side, and equation (D29) is achieved: 


P> 


16a 


(D29) 


Equation (D29) determines Stability Region 3. 

The graphs for the stability regions in figures Dl, D2, and D3 are limited to 0 < a < 2 andO < j3 <2 base on 
equations (Dll) and (D18). 



Figure D 1 . Stability Region 1 . 
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Figure D2. Stability Region 2. 



Figure D3. Stability Region 3. 

By superimposing all stability regions, one overlapped stability region is defined and shown in figure 7. Since 
many assumptions has been made during the development process, only conservative values for a and [> 
should be chosen. 
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APPENDIX E 
A Thermal Profile 

A thermal profile is very important in thermal tests. It consists of at least one segment that has a straight line 
relationship between temperature and time. In each segment, there are start time, end time, start temperature, 
and end temperature. The segment slope in °F/s is computed as follows: 

, end temperature - start temperature 

slope = 

end time - start time 

and the profile temperature is calculated with t is THRM PET in seconds: 

profile temperature = slope * (t - start temperature) + start temperature 
Table El below describes how a typical thermal profile is created for each zone. 


Table El. Thermal profile used in implementing GH closed-loop thermal control. 


Segment 

number 

Start time, 
second 

End time, 
second 

Start 

temperature, °F 

End 

temperature, °F 




Ambient 

80 

1 

0 

100 

80 

230 

2 

100 

200 

230 

530 

3 

200 

300 

530 

980 

4 

300 

400 

980 

1280 

5 

400 

500 

1280 

1480 

6 

500 

800 

1480 

1480 

7 

800 

1000 

1480 

1380 

8 

1000 

1200 

1380 

1180 
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